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Abstract 

In  this  paper,  we  consider  the  adaptive  detection  of  renewal 
processes  whose  inter-arrival  times  are  Gamma  distributed.  It 
is  shown  that  the  optimum  detector  exhibits  a two-dimensional 
estimator-correlator  structure  for  the  two  pertinent  parameters. 
When  the  underlying  statistics  are  partially  known,  the  estimates 
appearing  in  the  receiver  cannot  be  implemented.  Three  suboptimum 
schemes  with  surprisingly  good  small  sample  performance  are  derived 
and  compared. 
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. Introduction 

An  increasing  number  of  communication  ays terns  process  signals 
which  car.  be  modelled  as  point  processes.  These  occur  in  various 
aic-ic  such  as  optical  communications,  nuclear  medicine,  and  detec- 
tion of  seismic*  events.  Oftentimes,  the  sicnals  are  assumed  to  be 
roif.von  tirr.e-dt pendent  processes  and  detection  schemas  under  These 
assumptions  have  been  investigated  (fl]).  However,  many  processes 
depart  significantly  from  Poisson  statistics;  the  measure  of  depar- 
ture usually  is  taken  as  the  hazard  function  U2] ) which  is  con- 
stant under  the  Poisson  regime,  but  time-varying  fer  ether  renewel 
processes . 

A renewal  process  is  by  definition  a point  process  in  which 


the  sequence  of  times  between  occurence  of  events  consists  of  :.i.d. 
random  -ariables.  In  this  paper,  we  investigate  the  detection  c* 
renewal  processes  whose  inter-arrival  times  are  ?(u,k)  distributed, 
i . f . 

jf 

f(x!u,k)  - exp  (-  £ xjx1'*1^)  /r(k)  (I) 

With  two  parameters,  k and  a,  the  Gamma  distribution  is  a good 
r.adel  for  a variety  of  problems.  It  conveniently  describes  the 
* Poisson  regime  for  k * 1 and  measures  the  departure  from  Poisson 


statistics  through  the  parameters  k([2}).  In  particular,  charac- 
teristics of  bunching  are  quite  well  described  since 


!rom  (2),  it  cm  be  seen  that  if  k is  greater  than  on*,  vc  have 
spreading  of  the  observations  (i.e.  events  are  spaced  regularly 
•.round  the  mean  in  time)  whereas  if  k is  less  than  one,  events 
exhibit  a bunching,  or  correlated,  pattern. 

Ke  will  investigate  the  following  two  hypotheses  and  3^: 
...idt  r Hp,  noist  (dark  current)  is  received  and  the  process  is 
Poisson  with  mean  l/u0?  under  , the  observed  poin*-  process  con- 
tains a random  signal  to  be  detected  and  the  inter-arrival  times 
are  governed  by  (1).  We  will  assume  that  the  random  si-.jj.al  unde r 
!3^  modulates  the  information  bearing  parameters  u ar.d  V,  ?c  tV,:  t. 
theti  are  to  be  considered  as  random  variables.  Alter natively , 
one  right  consider  u as  the  information  bearing  parameter  while 
k reflects  the  unknown  dead  time  characteristic  of  a photomulti- 
plier device. 

In  order  to  determine  the  structure  of  the  optimum  detector 
minimizing  the  average  probability  of  error  under  a Bayesian 
criterion  or  maximizing  the  power  for  a fixed  probability  of  false 
alarm,  it  is  convenient  to  exploit  the  property  that  the  Gar.^a 
distribution  belongs  to  the  exponential  family.  Indeed,  let 


where 


and 


9 A (^.Bj)' 


0,  A-* 

1 u 


>2  * * 


h (x) 


Al 


(3) 


b (0  ) ^ logr(02)  - 02log(-01) 


Then,  (1)  car.  be  written  as 


f (x|0)  = h(x)exp(01x  + 02logx  - b (0  ; ) (4 

which  is  the  usual  exponential  form. 

In  the  next  section,  we  will  extend  some  of  thr*  rosilts  oi 
< f 3 1 ) to  the  two-dimensional  exponential  family  and  show  that, 
independent  cf  the  bivariate  prior  density  n(01,e2),  the  marginal 
density  f(x)  is  completely  determined  by  th*.-  conditional  mean 
estimates  (CMC)  of  0^  and  fi^.  This  resulting  form  for  the  margi- 
nal density  leads  to  a general  estimator-correla*  **r  structure  for 
detectors  based  on  likelihood  ratios. 

Since  the  optimum  detector  usually  cannot  be  implemented  be- 
cause o:  insufficient  a-priori  knowledge  cf  the  strtistics  of  u 
3nd  k,  we  will  investigate  the  properties  of  scrr.o  related  rub- 
optimum  detectors.  This  is  done  in  Section  III.  In  particular, 
we  will  utilire  a modified  and  a discrete  maximum  likelihood  esti' 
mate  ([4])  in  forminq  suboptimum  detectors.  The  simulations*  to 
be  discussed  illustrate  the  attractiveness  of  this  suboptimum  ap- 
proach, especially  in  the  important  small  sample  case. 
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'I,  Detection  of  a Renewal  Process  with  Gamma  Jnter-A r.  i.va  1 Times 
Bayesian  test 

We  suppose  that  under  both  hypotheses  ana  iij , n observa- 
tions (x . . i*l , . . . ,r. ) independent  and  identically  distributed 
fi.t.d.)  a.o  governed  by  (4?;  under  HQ,  6 = 9^  = f3,^,02u'*  is  3 
kno*»r.  vector,  whereas  under  Hj , £ is  a random  vector  v/it.n  bivariate 
>rijr  6 ~ r (6  t - r(6.,8^).  Moreover  we  assume  that  I!n  and  !?. 
o:eur  with  priors  equal  to  p0  and  p ^ respectively.  The  detect  io;i 
y.-cblen  adr.  its  a sufficient  statistic  ( f5] ) 


where1 


<W 


n 

E 

i=l 


(5) 


t 


2 


n 

- logx. 
i«l  1 


so  that  H0  and  Hj  become  equivalent  to  the  following.  Under  both 
hypotheses 

t ~ f'tjS)  - expffl^t^  + ®2t*'  ” nb  f® ) ♦Pit))  (€•) 

where  6 = 0^  is  a known  vector  under  HQ  and  under  I!j,  5 is  a rar.lr:. 
vector  with  prior  n (0  ) . Denoting  the  marginal  of  c ur  dvr  by 
fft),  the  optimum  detector  is  the  likelihood  ratio 


Lit) 


-fLU 

fTtfe 


?] 

7l 


£o 

Pi 


(7) 


f (t)  * J/f  (t | e)ti(e) de 


(8) 


I 


c 


Now,  f(t)  is  given  by 
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i.o. , 


f(t)  = exp(B(t))j*/exp(01t1  + 02t2  - nbr0)m(51,02'd01d52 

Take  the  partial  derivative  of  f(t)  with  respect  to  t^  ■■'.nd  t2: 

;j0.f(t|9)n(9)de 


Thus 


c loql 
fit 


mu  4 1 

hti 


l -•  1,2 


* E(9.|t)  * , . li2 


Since 


dlogf  (t ) 


dtx  + 


i 1 ' 

> 


i?on  substituting  (11)  into  (12),  one  obtains 

A A 

dlogf(t)  = 01(t)dt1  + 0 2 (t)dt0  * d3(t) 


( 13) 


(13)  is  a complete  differential.  Therefore  if  we  integrate  along 
a path  such  as  represented  in  Fig.  1,  we  get 

A t a t A 

f(t)  = K(0)exp(j  Qj  (u)du^  + J 92(u)du2  ■+  b(t)i  (1< 


where  K (6 ) is  the  normalizing  constant  and  tp  is  chosen  arbitrarily 


Fig.  1 


L-*t 

r($)  ft  } ft1(u)du1  + } d2(u)du2  - 0lotx  - e2£)t2  (15) 

t0  t0 

Arter  substitution  of  (6)  and  (2?)  into  (7),  one  can  writ*;  the 
2 i ke  1 i hood- r a t io  us 

L(t)  * K(§  )exp  (nb  (0  Q ) )exp(r  (ft  ) ) {1*-) 

following  (6],  the  constant  appearing  in  (16)  car.  be  written  in 
a more  convenient  way.  We  multiply  (16)  by  f < t_ j ^ q ) ?r.d  integrate 
with  respect  to  t.  Since  f(t)  integrates  to  1,  we  obtain 

[K(t  )exp(nb  (0  Q ) )]-1  « Eh  (exp(r(0))j  (17) 

n0 

• ebstituting  (17)  into  (16)  and  taking  logarithms,  the  log- likelihood 
ratio  becomes 

l(h  * r($)  - Eu  (expfr  (9  ) ) ) (16) 

H0 

ihis  is  compared  to  the  threshold  tn(pQ/p^)  for  an  optimum  Bayes 
test.  As  seen  from  (18)  and  (15),  the  optimum  receiver  i*  com- 
pletely determined  by  the  CME's  of  0^  and  ?2  and  displays  an  es*.’- 
r.itor-correlator  structure, 
o.  Neyman-Pearson  test 

l’he  Neyman-Pearson  test  is  easily  derivable  from  (18),  . 

A H.  A 

r (•  ) >X  v (0  ) (10) 

Ho 

A 

where  y (0 ) is  chosen  so  that  the  probability  of  false  alarm  is  set 
equal  to  a level  a. 

The  above  results  constitute  a canonical  detector  structure 
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'jr  renewal  processes  with  garura  intcr-arrival  ti.ni r.  tnc  r.:u.  t, 
of  course,  specify  the  prior  distribution  Whtn  thi3 

di  str -but ion  is  not  known,  Eqs.  (39'  and  f lr<  1 suggest  ror-Iacir.a  the 
CKE's  by  other  estimates  which  are  good  approximations  to  it  ::r.d 
which  require  less  prior  information.  Ihis  is  the  subject,  of  the 
r.axt  section. 

Finally,  it  should  be  clear  that  the  results  c 5 tl in  section 
os  easily  te  extended  to  the  n-dimensional  exponential  family. 


—“*«■■■■»  miSm, .innariim 
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IIT . Adaptive  Detection  of  Renewal  Processes 

Lot 

*10  ' - l/“0  • S20  * 1 
- - k/u  , e2  - k 


(20) 


Ylv-n , the  optimum  Bayesian  and  Neyraan-Ptarson  tests  ar-  given  by 
t lb  5 and  (10).  As  suggested  above,  these  tests  arc  otten  rot  used 
1 of  insufficient  prior  knowledge  or  beca.:.“e*  the  CM/'s  jf 

ar..i  S2  are  simply  difficult  to  implement.  Consequently,  it  is 
natural  to  investigate  the  properties  of  subopt imu^  detectors  c?b- 
* ained  by  substituting  suboptimum  estimators  for  th?  CME's  in  (l£  t 
•'c  (1°).  In  (|4J),  it  is  shown  that  good  approximation.,  to  CME's 
C3.:  h derived  from  modifications  of  the  maximum  likelihood  esti- 
r.ote  (MLE) . As  one  might  expect,  the  resulting  dctecscr  perfor- 
mance- is  close  to  the  optimum.  Vfhat  is  surprising  is  t.V  t fnis. 
is  true  even  for  very  small  samples  (n=3  or  4)  . .-‘e  new  dt-r  vc 

tx.rro  detection  schemes  based  on  the  MLE  and  modifications  of  it. 
This  is  done  in  increasing  order  of  assumed  prior  knowledge . The 
first  is  the  MLE  which  assumes  no  prior  knowledge  or.  p cr  k.  The 
truncated  MLE  assumes  that  the  range  of  p and  k is  know::.  Finally, 
the  discrete  MLE  assumes  further  that  the  parameter  k car.  only  take 
on  one  of  a finite  number  of  values. 

A.  MLE  Detector 

We  first  have  to  calculate  the  MLE 1 s of  p ana  k or,  equiva- 
lently, of  and  92.  From  (3)  and  (6)  the  maximum  likelihood 

equations  have  the  form 


tj/n 


(21) 


t(?2)  - log(-  9^)  * t2/r> 


T*  I I 7-v 


-10- 


where  u,  6^,  denote  the  MLE's  of  the  corresponding  parameters 
end  ♦ is  the  derivative  of  the  Gamma  function.  Th~  soluti .n  to 
(21)  is  not  immediate  and  does  not  lend  itself  to  analy- ic  ir.t  ’• 
oration.  However,  if  one  assumes  that  k(i.c.  8^)  in  .rv.ffjcier.tly 
'.?rce-  :o  that  Stirlino's  formula  ( [7J ) can  be  used,  have 

«(?2)  - lo,  ?2 

2 

ar.d 

h h. 

'*i  n 

1 tl  t'»  1 

*2  “ W (log  (^)  - ~)  (22 

For  later  use,  we  make  the  following  obeer^atior  : 

1)  and  ?2  car.  now  be  integrated. 

2)  k is  a reasonable  estimate  since  it  is  always  positive,  a 
preperty  which  stems  from  the  fact  that  the  arithmetic  rear.  s 
laryt-r  than  the  geometric  mean.  We  now  have  to  calculate  tne  in- 


tegrals 


. t 

I (8^  ) * j*  3^  (u)du^ 


I (^2)  ^ j*  (u^du2 
tn 


(24  ) 


where  the  integrations  should  be  performed  along  a convenient  pnth. 
In  Appendix  A,  these  integrations  are  carried  out,  the  final  result 


being: 


iffj)  ♦ - - f log  flog  (j^)  - 


Note  again  that  in  (25),  the  sign  of  the  argument  raises  no  problem 


since-  it  is  positive.  We  then  obtain  the  MLE  Bayesian  director 
by  substituting  (25)  into  (19).  We  also  have  to  .-nlculaf e the 
quantity 


K ' - e 


■ownte  (26)  as; 


H0exP  - h - 5 log  flog  h - , 


^1 

exp  [n  (log  — — - —*•)}  •:  t_ 

K ' sHr  r. — t~  n'2  ,xp>-  nleo  r~  * ^ <”> 

• tl09  _i  . c 

t.  t, 

ny  a theoren  due  to  Pitman  ([8].  page  217),  t,  and  loo  -=■ br> 

i ' n n 

independent.  This  property  permits  the  factcrir.c  of  the  expecta- 
tion in  (27)  and  since  t^  is  Gamma  distributed  under  r , we  have 

Vxpt- n 109  ^ + * nil  { <;r! 


t,  n-1 


dt, 


OD 

1 1 1 r . -1 
^ u " nn  0 1 * 


so  that 


Hence,  the  MLE  yields  an  undefined  Bayesian  detector,  r.  phenomenon 
encountered  for  ether  classes  of  problems  (f9j.  Fee.  V.r<,  and  flO] , 
Sec.  3.4  of  Chap.  2). 

In  contrast,  the  Neyman-Pearson  test  is  well  defined  and 
derived  by  substituting  (25)  into  (19).  Dividing  bv  n,  we  get 

un.-iicucii-li,  . i_  Ii . M’» . ,¥. 


'.T5-  . » m . , 


Truncated  MIE  Detector 


Since  the  optimum  detector  is  determined  by  CKE's,  one  mirht 
expect  that  by  modifying  the  MLE  for  some  given  partial  a -priori 
knowledge  of  the  parameters,  the  resulting  estimates  will  be 
closer  to  the  CV.E's,  and  the  associated  detector  will  ex’. .cit  ^ 


performance  which  is  closer  to  that  of  the  optimum.  Tru*  will 
indeed  be  the  case.  In  this  sub-section,  we  assume  prior  knew^ 


lrdgu  of  the  dynamic  range  of  u and  k,  i.e 


known 


We  consider  the  following  estimates 


are  gi ven  by 


In  Appendix  B,  we  calculate  the  integrals  I (0 1 ) and  T (?2 ) defined 
as  in  (23),  (24).  With  the  following  definitions 
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ff(x)] 


f (x ) if  f(x)  a 0 


0 if  f (x ) < 0 


c ix ) 


1 if  x a 0 


0 if  x < 0 


*1  * JT  ' y 


2 n 


a , log  _i . 


we  shw 


loa  — 
* n 


log^ 


-(KSj^)  + I (0  2 ) ) - - kj  (logy^-1 ) - k/ _ (-v,  ) + x c(y1-i; 

+ 7 Uog(2kj logy1)Jc(-v^)-(log{2ku)lcgi1]c(-vj)} 

k 

♦ k,  (logy^  - log4u)+  - ^lYi  “ ^u>  + + k4Y2  (32) 

- kt(y2"v*)+  “ 1 il°9^2k4an  x c^2^i-) 


+ j [log(2k^a)lc(y2-vu)  + ku  (y2  - vj 


Equation  (32)  is  obviously  a complicated  expression.  However,  it 
is  easily  imp lemen table  on  a computer  using  the  built-in  positive 
difference  function  or  on  special  purpose  hardware  u^ing  limiters. 
We  consider  only  the  associated  N-P  detector  which  is  obtained  by 
substituting  (32)  into  (19). 


*4V  - ’.v  r * 4 «•  - 

■ H 9*  1m  W ’ll?  j*  , ' 


the  a-priori  knowledge  of  * he  dyn:»  :c  ranae 


Here 


und  fisc  suppose-  that  k can  only  take  on  an  integer 


nw'» 


The  notation  here  should  not  be  confused  with  the  C>*. 
notation. 


The  resulting  N-P  detector  obtained  by  substituting  3 (6^)  and 

A 

I(e2)  into  (1*0,  is  derived  in  Appendix  C.  We  introduce  the 


following  notation: 


y2 


n 


vl  * lo«l  - 2 (X^i  )-l  ' * * 1 *u  ' \ 

"i  - e*pty2  + • 1 ■ 1 ‘ k> 


We  then  have 
I ) If  y.  > 1 


J (i  ♦ I(92))  * - kj'logyj-l)  + k,  (logy1  - Jogi-u)‘ 


k -k,  k.  -k, 

u l . u l . 

+ t (-v4)  + k^y2  + t (y,  - v.) 


i-1 


i-1 


2)  if  < yx  < l 


k -k. 
u l 


3)  If  yx  < 


~ (Kfl1)  + I(62)) 


k -k, 
u t 


(36 


£ (I(fl1)  + I(02))  - - kt(logy1-l)  (y2-v.)  + + kty2  131 


k -k, 

U ( B. 


- t (-A  - 


i-1 


- ** <iogut-D  + £ <y2-io?..t  ♦ 

k -k, 

+ 1 U \ X 

» S ("i-Vl' 

i-1  1 1 


(38 


Hi 


♦ tK  <“*  - yl>  + V2 


f «-f  ^ f i 


As  commented  on  in  Section  XII. B,  this  receiver  is  also  not  that 
difficult  to  implement. 

D.  Simulation  Results 

Simulations  have  been  performed  for  the  Nevman-Pearson  tests 
associated  with  V , 0 and  0,  and  are  denoted  respectively  by  DET.l. 
Dm. 2 and  DET.3.  Under  HQ,  the  observations  are  exponentially 
distributed  with  mean  l/uQ.  Under  • they  are  1'  • V./m  , k ) distri- 
buted, k is  uniformly  distributed  on  the  integers  Ck  ,k^+l, . . .k^} 
and  u is  independent  of  and  uniform  on  . For  this 

example,  the  optimum  test  which  we  designate  DET.4,  can  be  obtained 
d. rectly  from  the  likelihood  ratio  calculated  in  Appendix  D.  It 
should  be  noted  that  although  available  in  this  example,  this  de- 
tector cannot  be  set  into  an  estimator-correlator  structure  and, 
as  indicated  in  Table  3 below,  the  computing  time  required  for  its 
implementation  is  much  larger  than  that  of  any  of  the  tests  pre- 
viously described. 

We  simulated  hypotheses  HQ  and  1000  times  (m=l000)  and 
calculated  the  empirical  distributions  of  the  four  tests  under 
both  hypotheses.  To  determine  the  various  thresholds  for  a signi- 
ficance level  a,  we  used  the  following  non-parametric  method  dis- 
cussed by  Davis  and  Andreadakis  [11],  and  which  can  also  be  found 
in  ([12]).  Let  r(l),  r (2 ) , . . . ,r  (m)  be  the  order  statistics  of  any 
of  the  tests  investigated  under  HQ.  The  (l-a)^  quantile  qj  is 
such  that 

PrH0lr  * Ol-a1  - a 

Consider  the  event 

E - (r(j)  > q^} 


1 
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E occurs  if  at  least  (m-j+1)  values  of  r are  greater.  that  q 
corresponding  to  the  probability  of  having  at  least  (m-j+1) 


l~a ' 


successes  in  m Bernouilli  trials  with  a being  the  probability  of 


a success.  Hence 


Pr{E)  - Ift  (m-j+1. j) 


where  Ia(a,b)  is  the  incomplete  Beta  function.  In  this  case,  it 
can  be  approximated  by 


V vma  (1— c ) 


where 


a -- 

N (a,l)  « J e 2 dx 

J2h  -- 

Consequently,  for  m-1000,  if  we  choose  j»963  (or  j*918),  there  is 

a 96.4  4 probability  that  the  false  alarm  is  less  than  54  (or  1016) 

when  the  thresholds  are  taken  to  be  r(963)  and  r(9l8),  respectively. 

The  results,  summarized  in  the  following  tables,  illustrate 

some  significant  differences  in  the  small  sample  case  for  various 

values  of  the  parameters.  In  general,  DET.3  is  superior  to  DET.2 
which  in  turn,  performs  better  than  DET.l.  DET.3  is  quite  fre- 
quently much  better  than  DET.l  and  very  close  to  the  optimum.  In 
the  large  sample  case  (n  greater  than  10),  as  one  might  expect, 
the  detectors  have  similar  power. 


I 

1 


18 


Number 

of 

Samples 

M0 

W4 

Mu 

k4 

k 

u 

Power  of 
DET.l 

Power  of 
DET.2 

Power  of 
DET.3 

Power  of 
DET.4 
(optimum) 

3 

.3 

.5 

3 

1 

3 

.731 

.773 

.862 

.876 

4 

.5 

1 

3 

1 

3 

.689 

.784 

.818 

.658 

A 

*» 

3 

1 

5 

5 

9 

.354 

.443 

.52? 

.43  5 

4 

1 

.5 

3 

4 

7 

.433 

.484 

.615 

.681 

4 

1 

.5 

3 

2 

5 

.346 

.324 

. 4i 

.523 

10 

1 

.5 

3 

2 

5 

.708 

.797 

. 7C0 

.832 

Table  1.  a - 5% 


J 

) 

Number 

of 

Samples 

U0 

*u 

k< 

ku 

Power  of 
DET.l 

Power  of 
DET.2 

Power  of 
DET.3 

Power  of 
DET.4 
(optimum) 

d • 

i 

3 

.3 

.5 

3 

1 

3 

.816 

.824 

.911 

.915 

■ 

4 

.5 

1 

3 

1 

3 

.760 

.832 

.867 

.901 

• 

F 

M 

4 

3 

1 

5 

5 

9 

.549 

.699 

.726 

.760 

4 

1 

.5 

3 

4 

7 

.635 

.681 

.746 

.799 

|l. 

r 

ti 

4 

1 

.5 

3 

2 

5 

.495 

.501 

.5^7 

.660 

10 

1 

.5 

3 

2 

5 

.800 

.879 

.84? 

.893 

Table  2.  a ■ 10 X 


I 

I ; 


{ 


DET.l  DET.2  DET.3 

112  112  115 


DET.4 

495  seconds 


C 


Table  3. 


Approximate  Sum  of  Computing  Time 
for  the  First  and  Two  Last  Rows  of  Table  1. 


ML  ,,i4.  ttf-  . _*»  - . «V  • ' 'V  • 
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IV.  Conclusion 

In  this  paper,  we  investigated  the  detection  of  renewal  pro- 
cesses whose  inter-arrival  times  are  Gamma  distributed.  We  first 
developed  the  structure  of  the  optimum  Bayesian  and  Neyman -Pearson 
tests  for  the  two-dimensional  exponential  family.  The  main  charac- 
teristic of  these  detectors  is  that  they  fall  into  the  category  cf 
estimator-correlators  since  they  are  determined  by  integrals  of 
thr-  CME’s  of  the  two  pertinent  parameters.  One  implication  of  this 


structure  is  the  implementation  of  suboptimum  tests  by  substituting 


various  estimates  for  the  CME's. 

We  then  applied  the  estimator-correlator  property  to  the  case 
of  Gamma  distributed  observations  and  investigated  three  related 
tests.  The  first  detector,  DET.l,  is  based  on  the  NILE,  and  as  pre- 
viously observed,  the  Bayesian  test  is  undefined  whereas  the  Neyman- 
Pearson  version  seems  to  perform  quite  well  even  for  a small  number 
of  samples.  The  second  test,  DET.2,  is  based  on  the  truncated  MLE 
which  assumes  knowledge  of  the  dynamic  range  (boundary)  cf  the 


parameters.  Finally,  we  investigated  the  properties  of  0ET.3,  the 
test  based  on  the  discrete  MLE  of  k assuming  that  k can  only  assume 
a value  on  a finite  set  of  integers.  This  test  is  well-suited  for 
the  situation  where  observations  are  taken  at  the  output  of  a photo 
multiplier  with  a fixed  dead-time  characteristic.  DET.3  and  DET.2 
outperformed  DET.l  and  often  in  the  small  sample  case,  the  perfor- 
mance of  DET.3  is  markedly  superior  to  that  of  DET.l  and  very  close 
to  the  optimum.  Consequently,  based  on  these  preliminary  simula- 
tions, DFT.3  is  a fairly  adequate  test  for  the  detection  of  a large 
class  of  renewal  processes. 
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Since 


*1  *2 

1o9  ~i  - • n * 0 


the  admissible  points  -*  are  located  below  the  curve 

n 


y - log  — 


and  therefore  one  should  be  careful  in  choosing  the  path  of 
ir.ies  ration. 

t, 

1)  Case  1:  — 1 - > 1 

n 

We  integrate  along  path  (a).  Then 


1 (T, ) « J IT.  (u)du.  ■ - \ ^ 
1 A 1 * t 


U1  t2  dUl 
0 Ul  (109  — - -£) 


(3. A) 


Along  the  part  of  the  path  for  which  the  integral  does  not  vanish. 


we  have 


Thus 


*2  * <=20  " 0 


!l/n 


t P d?  _ n (•  ds 

1 (V  “ 2 J Hogo  " 2 { ^ 


I ) * - §“  l°9log 


(4. A) 


I ( 6*2 ) 


/ V,  (u)du,  - 


J u 

*0  log  — ^ 
n 


(5. A) 


I I 

I 


i -I 


Along  the  part  of  the  path  where  the  integral  does  not  vanish,  we 
have 


ui  p h 


Appendix  B 


Here,  we  integrate  the  estimates  5,  and  6 _ given  in  (311 


that 


ake  cn  is  located  below  the  curve  y * log  — . Several  rase 
be  investigated  which  will  be  referred  to  in  Fig.  B below 


Kc  assume  that  pu  > e,  k^  > ^ and  make  use  of  the  notation  intro- 
duced for  (32)  and  in  Fig.  B.  It  is  easily  verified  that 


5c  •*  k^  iff  y2  < 


k « k iff  * y2  * vu 


(I.e: 


K ■ k„  lff  *2  * vu 


- • C??o  1:  > 1 

1 

2k* 

a)  z,  * e < Yx  * bu 
Integrate  along  the  path  (a): 

yl  k 

Kff^/n  * - J ^ duj  - - k^  (logyj-1) 
e 1 


(2.B 


We  now  integrate  ?2  and  several  subcases  have  to  be  considered. 


i)  If  y2  > vu#  then 


Vt  . vu  du,  y2 

I (&2 )/n  * S k^du2  + j J logyi_U2  + J kudu2 

4* 


or 


I (»2)/n  - V«  * I 1o9  <i£>  * ku  0’2-vu>  <3' 


ii)  II<(‘J2<vu 


I (?2 )/n  * kjV^  - log  (2k^a) 


(4.: 


iii)  If  y2  < v^,  then 


I (ff2)/n  - k^y2 


(5J 


hen,  (3.B  - S.B)  can  be  rewritten  in  a single  formula,  i.e.. 


(6.B) 


I(02)/n  = k4y2  - ^<y2-vt)+  - 7 log(2kta)c(y2-v^) 
+ \ log(2kua)c<y2-Vu)  + ky  (y2  - vu>  + 


b)  Ij  y j > uu,  integrate  along  path  (b): 


or 


I (5^/n  - - JU  ^ dut  + J 1 - ^ d-n 
e 1 « 


I(e  )/n  = - k (logy  -1)  - ^ (y. ) 

v U m..  1 V. 


In  this  case,  I{§2>  is  again  given  by  (6.B). 

c)  lf  zu  s i’i  ‘ 2/  » then  integrate  alor.g  path  (c), 

U k , yl  da. 

I (§\)/n  - - f du.  - y f — ; 

1 e U1  1 2 z.  ullo<3’Ji 


or 


l(01)/n  « - kt(^~  - l)  - I iog(2k<tlogy1; 


Again,  I (S‘2 ) is  given  in  (6.B). 

d)  If  1 < y1  < zu,  we  integrate  along  path  (d),  so  tli.it 


io1)/n  * - vdr  ~ 1}  * i lo?  k1 " ku(l0?yi - dr} 

<■  u u 

and  again,  I(i»2)  is  given  by  (6.B).  We  can  include  (2.B,  7. 
9.B)  within  a single  formula,  i.e., 

Me^/n  - - k^dogyj-l)  - k^(-v^)+  + | log  (2k^logy1  )c  (-v^ ) 

- \ lon(2kulogy1)c(-vu)  + k (log  £i)+  - ^(y^-u  ) + 

u Mu  1 u 


(8.B) 

(9.H) 

B.  8 . B, 


(10.  B) 
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2.  Case  2:  y^  < 1 (Path  (e)) 
a)  If  y2  < vt 

Kff^/n  - - k^dogyj-1) 


b)  If  v4  s y2  s vu. 


I (ffj)/n 


- k^  (logsy -1 ) 


-hY 


du. 


2 s ui(1°9ui”y2> 
4 


KffjJ/n  * - k^dogs^-1)  - ^ log(2k^a) 


(11  .B) 


(12. B) 


c)  If  y2  > vu, 

l k. 

Kf^/n  « - k^dogs^-1)  - - log  - ku(vu»y2)  (13. B) 

Eqs.  (11. B - 13. B)  may  be  summarized  as 

I(B1)/n  - - k^dogyj-1)  - *4(y2“V+  - | log(2k,a)  (y2~vt)  + 

1 ♦ . (14.B) 

+ 2 log(2kua)  (y2-vu)  + ku (y2"vu) 

For  a),  b)f  c),  of  case  2,  we  have 


I (02  )/n  = k^y2 


(15.B) 


From  (6.B),  (10. B),  (14. B)  and  (15. B),  one  obtains  Eq.  (32). 


i 


_ 


r 


Appendix  C 


In  this  appendix,  the  discrete  MLE  detector  is  derived.  The 

AAA  A 

estimates  4,  k,  01  and  02  are  given  in  (33),  (34)  and  (35)  and 
have  various  forms  according  to  the  position  of  t in  the  plane. 
The  paths  of  integration  are  represented  in  Fig.  Cl  and  Fig.  C2 
below. 


logy 


Fig.  Cl 


..  r»r-  -3T  -•  ' Lg T ~ 7>P  Ar  * , 


v;«  assume  that 


< 1 


U > 

u 


and  i 1 


With  the  notation  introduced  below  (35),  it  is  readily  verified 
that  (34)  is  equivalent  to: 

ft  = k. 


iff  y2  * Vj 


ft  * k^+i  iff  < y 
ft  « k 


2 * Vi+1  ' i 


i,...,k,  - v - i a.c) 

a <• 


u 


iff  < y. 


2.  Case  1 : If  > 1,  then  consider  Fig.  c.l  and  integrate  along 
the  appropriate  path. 

A 

a)  If  yj^  a z^,  1(0^^)  is  given  in  (2.B)  and  (7.B). 

i)  If  y,  s v., 

fc2 

I(§2)  - J ktdu2  * nk*y2  U.O 


ii)  If  vx  < y2  * v2, 

I (§2 )/n  * k^Vj  + (k^+1)  (y2-v1)  - k^y2  + (y2“'»1)  <3.0 

iii)  If  Vj  < y2  < vJ( 

I(ft2)/n  « k^y2  + (y2~v ^ ) ♦ (y2-v2)  (4.0 

Cne  can  summarize  (2.C  - 4.C)  and  the  other  subcases  as 

k -k. 

A Ml. 

I(fi2)/n  - k^y2  + 2.  (y2"vi)  (5.0 


b)  If  *2  < Yi  < *,i  then 


I : 


) is  still  given  by  (5.C) 


Cases  a),  b),  c) . and  all  ether  subsequent  cases  for 


I (8.  )/n  = - k.  (logy.-l)  + k.  (leg 


while  1(8.)  is  given  by  (5.C) 


< y.  < 1,  then  consider  Fig.  C.2  and  integrate 


along  the  appropriate  paths 


fc.  (logy.-l)  + 
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•?.  case  3 : If  y 1 s p^,  we  have  to  consider  the  various  valuer 


3. 

Case  3: 

tr.at  p^  can  < 

a) 

If  > : 

i) 

if  y2  > 

ii) 

If  s 2 < 

A 

1 (B 1 

(12 

.C),  (13. 

A 

I(01)/n 

b) 

If  s2  < 

I(e1)  /n  * - k^Uogp^-l)  + “ (M^-y^ 


V 


(12.C) 


IfB^/n  = - k/(loga/-l)  + jj-  (s^y^  ♦ ^ (13. C) 

(■  4 


k -k, 
u K 


IfB^/n  « - k^(logu^-l)  + j~-  *si“yl*+  + ^l4-c) 


If  s2  < y^»  we  have  using  (10. C) 

KB^/n  - - k4(logpt-l)  ♦ (y2-logut  + jkT+l0  - (y^,) 

v ^ 


k.  +1 


(15. C) 


and  in  general  for  ^ yl  C *i#  * * 2, . . . ,ku~k^ » we 

I(61)/n  - - kt(logut-l)  -t  (y2  - logp^  + j^) 
1 . k.+l 

+ u c (8i-yi)  r^~(Mryi) 

Ml  i«2  1 1 4 1 


(16. C) 


We  can  regroup  (14. C)  and  (16. C)  as 

. ♦ 

K^J/n  » - k^logp^-l)  + (y2  - logp^  + mT+I^ 

1 ♦ 1 V\  4 ^ *, 

- — - (s.-p. ) + ~-  E (s  ,-y. ) ♦ r^u,- y.) 

ut  1 4 ut  icl  i 1 mj . * 1 


(17.C) 


■ 


•r%* 


*1*7.  «£V 
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Ir.  fact,  (17. C)  can  be  generalized  to  the  cases 
si+^  < * s^  , i * 2,...,ku»k^  as  follows: 

k -k. 

A U t 

Kfl^/n  « - k^  (logu^-1)  + £ (y,  - logu,  + 


l 2 


)' 

(k^+i  )-l 


-57  [ “^(W+-  ,si-yi,+] + ^-yi 


(18. C) 


For  < 1,  we  always  have 


I (80  - nk.y. 


(19. C ' 


Finally,  Eqs.  (5.C),  (8.C),  (11. C),  (18. C)  and  (19. C)  yield  the 
results  stated  in  (36),  (37),  and  (38). 
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Appendix  D 


We  derive  the  optimum  detector  DET.4.  The  marqinal  of  t 
under  Hj  which  we  denote  by  f(t),  can  be  written  as: 


f(t)  « JJf  (t|u,k)n(k)n(u)dkdu 


(1.3) 


where  ir(k)  and  tt  (u ) are  the  priors  on  k and  u respectively.  Con- 
sequently, we  have 

<tt’"  i:  H>nk f4txp'-  - ws,t,,4‘ 


f2.D) 


J - Ju  «<p<-  l tt)  -if  4. 


(3.D) 


We  integrate  by  parts: 


Iterating  the  integrations  by  parts  and  denoting 


14. D) 


Hence,  substituting  (5.D)  into  (2.D)  and  using  (6),  the  likelihcod- 
rc.tio  which  we  denote  by  L(t),  is  equal  to* 


L(t)  4 gftfe  ) “ c 9X^T  ~ n'k'-l!n exp[  (v-l)t  ) 

1 1 <r  u0  k«k^  tj  Tft)  2 


-z  nk"2  z1 

e r fr 

i-0  11 


kt. 


u 


kt. 


(6.3) 


where 


u. 


n 


c ’ W' 


-35- 


References 


1.  Snyder,  D.,  Random  Point  Processes,  Wiley,  N.Y..  1975. 

2.  Cox,  0.  and  Lewis,  P.,  Tha  Statistical  Analysis  of  Events. 
Methuen,  London,  1966. 

3.  Schwartz,  S.C.*  "Conditional  Mean  Estimates  and  Bayesian 
Hypothesis  Testing,"  IEEE  Trans.  Inform.  Theory,  Vol.  IT-21, 
Nov.  1975,  pp. 663-665. 

4.  Fogel , A.  and  Schwartz,  S.C.,  "On  Approximations  to  Condi- 
tional Mean  Estimates,"  Proc.  Fourteenth  Annual  Allercon 
Conference  on  Circuits  and  Systems  Theory,  Oct.  1976, 

pp. 944-953. 

3.  Ferguson,  T.,  Mathematical  Statistics,  Academic  Press,  1967. 

6.  Schwartz,  S.C.,  "The  Estimator-Correlator  for  Discrete-Time 
Problems,"  IEEE  Trans.  Inform.  Theory.  Vol.  IT-23,  Jan.  1977, 
pp. 93-100. 

7.  Abramowitz,  M.  and  Stegun,  I.,  Handbook  of  Mathematical 
Functions,  Dover,  N.Y. 

8.  Pitman,  E.,  "The  Closest  Estimates  of  Statistical  Parameters, 
Proc.  Camb.  Phil.  Soc..  Vol.  33,  1937,  pp. 212-222. 

9.  Fogel,  A.  and  Schwartz,  S.C.,  "Bounds  on  Probability  of  Error 
and  Suboptimal  Detectors,"  submitted  to  IEEE  Trans.  Inform. 
Theory. 

10.  Fogel,  A.,  "Adaptive  Detection  and  Estimation,"  Ph.D.  Thesis, 
Princeton  University,  Feb.  1978. 

11.  Davis,  M.H.A.  and  Andreadakis,  E.,  "Exact  and  Approximate 
Filtering  in  Signal  Detection:  An  Example,"  IEEE  Trans . 
Inform.  Theory,  Vol.  IT-23,  Nov.  1977,  pp. 768-773. 

12.  Kendall,  M.  and  Stuart,  A.,  The  Advanced  Theory  of  Statistics 
Vol.  2,  Griffin,  1961. 


v ***-  -V* 


-36- 


OFFICE  OF  NAVAL  RESEARCH 
STATISTICS  AND  PROBABILITY  PROGRAM 

BASIC  DISTRIBUTION  LIST 
UNCLASSIFIED  TECHNICAL  REPORTS 


JULY 

Copies 

Statistics  and  Probability 
Program  (Code  435) 

Office  of  Naval  Research 
Arlington,  VA  22217  3 

Director,  Naval  Research 
Laboratory 

Attn:  Library,  Code  2029 
(GNRL) 

Washington,  D.C.  20375  6 

Defense  Documentation  Center 
Cameron  Station 

Alexandria,  VA  22314  12 

Technical  Information  Division 
Naval  Research  LaDoratory 
Washington,  D.C.  20375  6 

Office  of  Naval  Research 

New  Ycrx  Area  Office 

715  Broadway  - 5th  Floor 

Attn:  Dr.  Robert  Grafton 

New  York,  New  York  10003  1 

Director 

Office  of  Naval  Research 
Branch  Office 

666  Summer  Street , Bldg.  114,  Sec.  D 

Attn:  Dr.  A.L.  Powell 

Boston,  NA  02210  1 

Director 

Office  of  Naval  Research 
Branch  Office 
536  South  Clark  Street 
Attn:  Deputy  and  Chief  Scientist 
Chicago,  Illinois  60605 


1977 

Copies 

Director,  Office  of  Naval 
Research  Branch  Office 
1030  East  Green  Street 
Attn:  Dr.  A.R.  Laufer 
Pasadena,  CA  91101  1 

Director,  Office  of  Naval 
Research  Branch  Office 
1030  East  Green  Street 
Attn:  Dr.  Richard  Lau 
Pasadena,  CA  91101  1 

Office  of  Naval  Research 
San  Francisco  Area  Office 
One  Hallidie  Plaza  - Suite  601 
San  Francisco,  CA  94102  1 

Office  of  Naval  Research 
Scientific  Liaison  Group 
Attn:  Dr.  Bruce  J.  McDonald 
Scientific  Director 
American  Embassy  - Tokyo 
APO  San  Francisco  96503  1 

Technical  Library 

Naval  Ordnance  Station 

Indian  Head,  Maryland  20640  1 

Naval  Ship  Engineering  Center 
Philadelphia 

Division  Technical  Library 
Philadelphia,  PA  19112  1 

Bureau  of  Naval  Personnel 
Department  of  the  Navy 
Technical  Library 
Washington,  D.C.  20370  1 
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Library,  Code  0212 

Naval  Postgraduate  School 

Monterey,  California  93940  1 

Library 

Nava',  Ocean  Systems  Center 

San  Diego,  California  92152  1 

App'.ied  Mathematics  Laboratory 
David  Taylor  Naval  Snip  Research 
and  Development  Center 
Attn:  Mr.  Gene  K.  Gleissr.er 
Bethesda,  Maryland  20084  1 

Stanford  University 
Department  of  Operations  Research 
Attn:  Prof.  G.  Lieberman 
Stanford,  CA  94305  1 

Florida  State  University 
Department  of  Statistics 
Attn:  Prof.  R.  Serfling 
Tallahassee,  FL  32306  1 


Florida  State  University 
Department  of  Statistics 
Attn:  Pi  of.  R.A.  Bradley 
Tallahassee,  FL  32306  1 

Princeton  University 
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Columbia  University 
Department  of  Civil  Engineering 
and  Engineering  Mechanics 
Attn:  Prof.  C.  Derman 
New  York,  New  York  10027  1 

Columbia  University 

Department  of  Mathematics 

Attn:  Prof.  H.  Robbins 

New  York,  New  York  1 

New  York  University 
Courant  Institute  of 
Mathematical  Sciences 
Attn:  Prof.  W.M.  Hirsch 
New  York,  New  York  10453  1 

University  of  North  Carolina 

Department  of  Statistics 
Attn:  Prof.  W.L.  Smith 

Chapel  Hill,  NC  27514  1 

University  of  North  Carolina 
Department  of  Statistics 
Attn:  Prof.  M.  Leadbetter 

Chapel  Hill,  NC  27514  1 

University  of  California 
San  Diego 

Department  of  Mathematics 
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Department  of  Statistics 

Attn:  Prof.  O.W.  Tukey 

Princeton,  NO  C3540  1 

Pr i r.cston  Ur. i vers i ty 
Department  of  Statistics 
Prof.  G.S.  Watsor, 

Princeton,  NJ  08540  1 


Attn:  Prof.  M.  Rosenblatt 
La  Jolla,  CA  92093 

State  University  of  New  York 
at  Buffalo 

Statistical  Science  Division 
Attr.:  Prof.  E.  Parzen 
Amherst,  New  York  14226 


Stanford  Un i ve  r s i ty 


University  of  California 
Operations  Research  Center 
College  of  Engineering 
Attn:  Prof.  R.E.  Barlow 
Berkeley,  CA  94720 

Yale  University 
Department  of  Statistics 
Attn:  Prof.  F.J.  Anscombe 
New  Haven,  CT  06520 
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Purdue  University 
Department  of  Statistics 
Attn:  Prof.  S.S.  Gupta 
Lafayette,  Indiana  47907 

Cornell  University 
Department  of  Operations  Research 
Attn:  Prof.  R.E.  Bechhofer 
Ithaca,  New  York  14850 

Southern  Methodist  University 
Department  of  Statistics 
Attn:  Prof.  D.B.  Owen 
Dallas,  Texas  75275 

Daniel  H.  Wagner,  Associates 
Station  Square  One 
Paoli,  Pennsylvania  19301 

Stanford  University 
Department  of  Operations  Research 
Attn:  Prof.  A.F.  Veinott 
Stanford,  CA  94305 

Stanford  University 
Department  of  Operations  Research 
Attn:  Prof.  D.L.  Iglehart 
Stanford,  CA  94305 

Stanford  University 
Department  of  Statistics 
Attn:  Prof.  Solomon 
Stanford,  CA  94305 

University  of  North  Carolina 
Department  of  Statistics 
Attn:  Prof.  C.R.  Baker 
Chapel  Hill,  NC  27514 

Clemson  University 
Department  of  Math.  Sciences 
Attn:  Prof.  K.T.  Wallenius 
Clemson,  SC  29631 

Case  Western  Reserve  University 
Department  of  Mathematics  and  • 
Statistics 

Attn:  Prof.  S.  Zacks 
Cleveland,  Ohio  44106 


Yale  University 
Department  of  Statistics 
Attn:  Prof.  I.R.  Savace 
New  Haven,  CT  06520 

Southern  Methodist  University 
Department  of  Statistics 
Attn:  Prof.  W.R.  Schucany 
Dallas,  Texas  75275 

University  of  Missouri 
Department  of  Statistics 
Attn:  Prof.  W.A.  Thompson,  Jr. 
Columbia,  Missouri  65201 
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Rice  University 

Department  of  Mathematical  Sciences 
Attn:  Prof.  J.R.  Thompson 
Houston,  Texas  77001  1 

Naval  Postgraduate  School 
Department  of  Operations  Research 
and  Admini strati ve  Sciences 
Attn:  Prof.  P.A.W.  Lewis 
Monterey,  California  93950  1 

University  of  California 
Department  of  Applied  Physics 
and  Information  Science 
Attn:  Prof.  E.  Masry 
La  Jolla,  CA  92093  1 

University  of  California 

School  of  Engineering 

Attn:  Prof.  N.J.  Bershad 

Irvine,  California  92664  1 

University  of  California 
School  of  Engineering  and 
Applied  Science 
Attn:  Prof.  I.  Rubin 
Los  Angeles,  CA  90024  1 

Virginia  Polytechnic  Institute 
and  State  University 
Department  of  Industrial  Engineering 
and  Operations  Research 
Attn:  Prof.  R.L.  Disney 
Blacksburg,  VA  24061  1 
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Copies 

..dvc'i  Postgraduate  School 
Dc-p*.  r trier,  t of  Operations  Research 
arc  Ac...; .-.1  strati ve  Sciences 
Attn:  Prof.  J.D.  £sary 
.\onterey,  California  93940  1 

PoV/tec .*.r.\c.  Institute  of  New  York 
^ w . ov  a.ectric&i  w..gir.eering 
Attn:  Prof.  M.L.  Sr.oo.tan 
Brooklyn*  .New  York  11201  1 

Carnegie-Mel  lor.  University 
Department  of  Statistics 
Attn:  Prof.  J.B.  Kudsne 
Pittsburgh,  ?A  15213  1 

Coicracio  State  University 
Dept,  of  Electrical  Engineering 
Attn:  Prof.  L.L.  Sciiarf,  Jr. 

Fort  Collins,  CO  83521  1 

Virginia  Polytechnic  Institute 
anc  State  University 
Department  of  Economics 
Attn:  Prct . X.o.  H'.nich 
Blacksourg,  Virginia  24061  1 

Rockwell  International  Corporation 

Science  Center 

P.0.  Box  1035 

nttn:  - r . x.  R.  Mann 

Thoasar.c  Oaxs*  CA  9.3oC  1 

University  of  Southern  California 
Department  of  Quantitative 
Business  Analysis 
Attr.:  Prof.  W.R.  Blischke 
Los  Angeles,  CA  50007  1 

Southern  Methodist  University 
Computer  Sc*.ence/Gperations 
Research  Center 
Attn:  Prof  U.N.  3hat 
Dallas,  Texas  75275  1 

University  of  Florida 
Department  of  Industrial  and 
Systems  Engineering 
Attn:  Prof.  R.S.  Leavenworth 
Gainesville,  FL  32611  1 
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Massachusetts  Institute  of 
Technology 

Department  of  Mathematics 

Attr.:  Prof.  H.  Chernoff 

Cambridge,  MA  02139  1 

Desmat'ics,  Inc. 

P.0.  Box  618 

Attn:  Dr.  D.E.  Smith 

State  College,  PA  16801  1 

Cornell  University 
Department  of  Operation 
Attn:  Prof.  J.A.  Muckstadt 

Ithaca,  New  York  19850  1 

Commandant  of  the  Marine  Corps 
(Code  AX) 

Attn:  Dr.  A.L.  Slafkosky 
Scientific  Advisor 

Washington,  O.C.  20380  1 

Program  in  Logistics 

The  George  Washington  University 

Attn:  Dr.  W.H.  Marlow 

707  22 r.d  Street,  N.W. 

Washington,  D.C.  20037  1 

Navy  Library 

National  Space  Technology  Laboratory 

Attn:  Navy  Librarian 

Bay  St.  Louis,  Mississippi  39522  1 

National  Security  Agency 
Attn:  Dr.  James  R.  Maar 
Fort  Meade,  Maryland  20755  1 

Naval  Coastal  Systems  Laboratory 
Code  P761 

Attn:  Mr.  C.M.  Bennett 

Panama  City,  FL  32401  1 

National  Security  Agency 

Attn:  Mr.  Glenn  F.  Stahly 

Fort  Meade,  Maryland  20755  1 

Mr.  F.  R.  Del  Priori 
Code  224 

Operational  Test  and  Evaluation 
Force  (OPTEVFOR) 

Norfolk,  Virginia  23511  1 


V **  -f*  • v X—  * % r ■ r 

% , . . f ...  t . 


Copies 


Copies 


.laval  Electronic  Systems  Command 
(NELEX  320) 

National  Center  No.  1 

Arlington,  Virginia  20360  1 

Dr.  A.  Petrasovlts 

Room  2078,  Food  & Drug  Building 

Tunney’s  Pasture 

Ottowa,  Ontario  K1A-0L2,  Canada  1 

J.S.  Lee  Associates,  Inc. 

Attn:  Dr.  J.S.  Lee 

2001  Jefferson  Davis  Highway 

Suite  802 

Arlington,  VA  22202  1 

Arthur  D.  Little,  Inc. 

Attn:  Mr.  Wendell  6.  Sykes 
Acorn  Park 

Cambridge,  MA  02140  1 

Claremont  Men's  College 
The  Institute  of  Decision  Science 
for  Business  and  Public  Policy 
Attn:  Dr.  Jane*.  M.  Myhre 
Claremont,  CA  91711  1 

Cornell  University 
Sibley  School  of  Mechanical  and 
Aerospace  Engineering 
Attn:  Prof.  S.L.  Phoenix 
Ithaca,  New  York  14850  1 

University  of  Missouri 
Department  of  Statistics 
Attn:  Prof.  R.W.  Madsen 
Columbia,  Missouri  65201  1 

State  University  of  New  York, 
Binghamton 

Department  of  Math.  Sciences 
Attn:  Prof.  D.L.  Hanson 
Binghamton,  New  York  13901  1 

Virginia  Polytechnic  Institute 
and  State  University 
Department  of  Statistics 
Attn:  Prof.  L.D.  Lee 
Blacksburg,  VA  24061  1 


University  of  Pittsburgh 
Department  of  Mathematics 
Attn:  Prof.  H.  W.  Block 
Pittsburgh.  PA  15260  1 

University  of  North  Carolina 
Department  of  Mathematics 
Attn:  Prof.  M.  Abdel -Hameed 
Charlotte,  NC  28223  1 

Illinois  Institute  of  Technology 
Department  of  Industrial  and 
Systems  Engineering 
Attn:  Prof.  M.  Zia  Hassan 
Chicago,  Illinois  60616  1 

Kansas  State  University 
Department  of  Industrial  Engineering 
Attn:  Prof.  F.A.  Tillman 
Manhattan,  Kansas  66506  1 

Naval  Postgraduate  School 
Dept,  of  Operations  Research  and 
Administrative  Sciences 
Attn:  Prof.  D.  P.  Gaver 
M ortterey,  California  93940  1 

Lehigh  University 
Dept,  of  Industrial  Engineering 
Attn:  Prof.  G.E.  Whitehouse 
Bethlehem,  PA  18015  1 

Stanford  University 
Stanford  Electronics  Laboratories 
Attn:  Prof.  Bernard  Widrow 
Stanford,  California  94305  1 

Naval  Research  Laboratory 
Advanced  Projects  Group 
(Code  8103) 

Attn:  Mr.  David  A.  Swick 
Washington,  D.C.  20375  1 

Defense  Communications  Agency 
Defense  Communications  Engineering 
Center 

Attn:  Dr.  M.  J.  Fischer 
1860  Wiehle  Avenue 

Reston,  Virginia  22090  1 
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J.S.  Army  Research  Office 

P.0.  Box  12211 

Attn:  Dr.  J.  Chandra 

Research  Triangle  Park,  NC  27706  1 

Naval  Tea  Systems  Command 
(XSEA  03F) 

Attn:  Miis  B.S.  Orleans 
Crystal  Plaza  #6 

Arlington,  Virginia  20360  1 

Office  of  The  Director 

Bureau  of  The  Census 

Attn:  Mr.  H.  Nisselson 

Rm.  2025,  Federal  Bldg.  3 

Washington,  D.C.  20233  1 

Defense  Logistics  Studies 
Information  Exchange 
Army  Logistics  Management  Center 
Attn:  Mr.  J.  Dowling 
Fort  Lee,  Virginia  20390  1 

OASD  (I&L),  Pentagon 

Attn:  Mr.  Charles  S.  Smith 

Washington,  D.C.  20301  1 

The  American  University 
Dept,  of  Mathematics,  Statistics 
and  Computer  Science 
Attn:  Prof.  Joseph  Blum 
Washington,  D.C.  20016  1 

The  George  Washington  University 
Department  of  Operations  Research 
Attn:  Prof.  N.D.  Singpurwalla 

Washington,  D.C.  20052  1 

Stanford  University 
Department  of  Statistics 
Attn:  Prof.  D.  0.  Slegmund 
Stanford,  California  94305  1 

San  Diego  State  University 
Dept,  of  Electrical  Engineering 
Attn:  Prof.  F.J.  Harris 
San  Diego,  California  92182 


Indiana  University  Foundation 
Department  of  Mathematics 
P.0.  Box  F 

Attn:  Prof.  M.L.  Purl 
Bloomington,  Indiana  47401  1 

Polytechnic  Institute  of  New  York 
Dept,  of  Operations  Research  and 
System  Analysis 
Attn:  Prof.  L.H.  Herbach 
Brooklyn,  New  York  11201  1 

University  of  California 

College  of  Engineering 

Attn:  Prof.  S.  M.  Ross 

Berkeley,  California  94720  1 

Union  College 

Institute  of  Administration  and 
Management 

Attn:  Prof.  L.A.  Aroian 
Schenectady,  New  York  12308  1 

Southern  Methodist  University 
Department  of  Statistics 
Attn:  J.E.  Boyer,  Jr. 

Dallas,  Texas  75275  1 

University  of  Wisconsin 
Department  of  Statistics 
Attn:  Prof.  Grace  Wahba 
Madison,  Wisconsin  53706  1 

Naval  Underwater  Systems  Center 
Attn:  Dr.  Derrill  J.  Bordelon 
Code  21 

Newport,  Rhode  Island  02840  1 
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